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Abstract 

We present a systematic exploration of dark energy and modified gravity models contain¬ 
ing a single scalar field non-minimally coupled to the metric. Even though the parameter 
space is large, by exploiting an effective field theory (EFT) formulation and by imposing 
simple physical constraints such as stability conditions and (sub-)luminal propagation of 
perturbations, we arrive at a number of generic predictions. (1) The linear growth rate of 
matter density fluctuations is generally suppressed compared to ACDM at intermediate red- 
shifts (0.5 < z < 1), despite the introduction of an attractive long-range scalar force. This 
is due to the fact that, in self-accelerating models, the background gravitational coupling 
weakens at intermediate redshifts, over-compensating the effect of the attractive scalar force. 
( 2 ) At higher redshifts, the opposite happens; we identify a period of super-growth when the 
linear growth rate is larger than that predicted by ACDM. (3) The gravitational slip param¬ 
eter 77 —the ratio of the space part of the metric perturbation to the time part—is bounded 
from above. For Brans-Dicke-type theories p is at most unity. For more general theories, 
77 can exceed unity at intermediate redshifts, but not more than about 1.5 if, at the same 
time, the linear growth rate is to be compatible with current observational constraints. We 
caution against phenomenological parametrization of data that do not correspond to predic¬ 
tions from viable physical theories. We advocate the EFT approach as a way to constrain 
new physics from future large-scale-structure data. 
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1 Introduction 


Understanding the nature of the present cosmic acceleration is an important and fascinating 
challenge. The standard paradigm—a cosmological constant + Cold Dark Matter within the 
framework of general relativity (ACDM)—has so far held up remarkably well when tested against 
cosmological data 130' This is especially true for data on the background expansion history. 
Large scale structure data, in other words data that concern the fluctuations, are improving 
in precision—current constraints are broadly consistent with ACDM, although mild tensions 
exist (see e.g. [5j-{8] ). New surveys, such as DES j9], Euclid |10|, DESI 11 , LSST 12 and 
WFIRST fl3|, are expected to greatly tighten these constraints. 


13 


An important question is 


whether or how these observations can be used to distinguish ACDM from other dark energy 
models [I4|jl7| . 

From the theoretical perspective, any form of dark energy that is not the cosmological 
constant would have fluctuations. These dark energy fluctuations can couple to matter or not. 
Or, in the frame where matter couples only to the metric (the frame we use in this paper), 
these dark energy fluctuations can kinetically mix with the metric fluctuations or not. Models 
that have no such mixing are quintessence models—they generally give predictions for structure 
growth fairly close to ACDM, especially if the background expansion history is chosen to match 
to the observed one, with the equation of state index close to —1. Models that have such mixing 
are modified gravity models, our primary interest in this paper. 

Effective field theory provides a framework for systematically writing down the action that 
governs the dynamics of such dark energy fluctuations jl8]-27] (see 28 - 301 for a numerical imple¬ 
mentation of this formalism). The background cosmic expansion is treated as a given—matching 
that of the best-fit ACDM for instance—while the dark energy fluctuations are encoded by a 
single scalar, the Goldstone boson associated with spontaneously broken time-diffeomorphism 
in the gravity/dark energy sector. In this paper, we are interested in the linear evolution of 
fluctuations, thus we retain terms in the action up to quadratic order in a gradient expansion. 
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An important feature of the effective field theory for dark energy is that the Ricci scalar comes 
with a general time-dependent coupling. This allows the accelerating background expansion to 
be driven by modified gravity effects, as opposed to what resembles vacuum energy. 

At the level of linear perturbations, non-standard gravitational scenarios effectively result in 
a time- and scale-dependent modifications of the Newton’s constant G e fr and of the gravitational 
slip parameter r/ 1311. The former quantity captures information about the way mass fluctua¬ 
tions interact in the universe, while the gravitational slip parameter quantifies any nonstandard 
relation between the Newtonian potential $ (time-time part of the metric fluctuations) and the 
curvature potential T (space-space part). 

Models of dark energy contain in principle so many parameters that it might seem hopeless 
to come up with robust, generic predictions for G e ff, r/. For example, Horndeski theory 


32,33 


the most general theory containing one additional scalar degree of freedom </>, with second order 
equations of motion—depends on four arbitrary functions of 4> and of the kinetic term (<9</>) * 2 . If 
we relax the condition of no-higher derivatives in the equations of motion to that of absence of 
pathological ghost instabilities, we end up with an even larger beyond Horndeski set of theories, 
dependent on six arbitrary functions of <f> and {d<f>) 2 (34[[35] . 

One of the main purposes of this paper is to show that, despite such apparent freedom, the 
behavior of G e s(z) and rj(z) as functions of the redshiffQ has definite features common to all 
healthy dark energy models (models with no ghost and gradient instabilities, nor superluminal 
propagation) within the vast Horndeski class. As a corollary, we show that popular, phenomeno¬ 
logical parameterizations of G e fj (z) and r](z) do not capture their redshift-dependence in actual, 
physical models. 

The phenomenology of theories containing up to two derivatives in the equations of motion 
has been explored in Ref. 24 , hereafter paper I, where a complete separation between back¬ 


ground and perturbations quantities has been obtained (see also |37|) and a systematic study 
of the growth index of matter density fluctuations has been initiated. As shown in paper I, the 
background evolution and the linear cosmological perturbations for any scalar-tensor theory of 
this type are entirely captured by one constant and five functions of timej^] 


{^rn, H (t), m(*) 5 M 2 (t), M3 (t), e 4 (t)} 


( 1 ) 


The above ingredients allow one to span the entire set of Horndeski theories, which can be 
seen as generalizations of galileon-theories |40|. Since also models of massive gravity 1411 and 
bi-gravity 42 reduce to the galileons in the relevant (“decoupling”) limit [43,44], aspects of 


these models are captured by our analysis. In the above, = p m (to) / (3m,p 2 Hq) is the present 
fractional energy density of non-relativistic matter; H(t), the Hubble parameter, encapsulates 
the background expansion history; the fi s are parameters of the perturbation sector with the 
dimension of mass, typically of order Hubble: p(t) is the non-minimal coupling of Brans-Dicke 
( BD ) theories, M 3 appears in cubic galileon- and Horndeski-3 (. H3 ) theories, M 2 affects the sound 
speed of the scalar fluctuations, but has otherwise no bearing on the linear growth of matter 
fluctuations; e 4 is a dimensionless order one parameter present in galileon-Horndeski 4 and 5 
( HJf.5 ) Lagrangians. Note that the full Horndeski theory ( Hf5 ) includes H3 which, in turn, 


J The scale dependence of G e g and y arises from (possibly time-dependent) mass terms for the dark energy 
fluctuations. A natural value for the mass is the Hubble scale, implying essentially no scale dependence in G e g 
and y for fluctuations on scales much smaller than the Hubble radius. An observable scale dependence for the 
growth of such fluctuations can only arise if one introduces a mass scale higher than Hubble. This is the case for 
chameleon models for instance 36 . We will not consider this possibility here. 

2 With respect to the equivalent “a-parameterization” introduced in Ref. [37] (see App. 0, the “y- 
parameterization” used here is more theory-oriented: as summarized in the text, our couplings are in direct 
correspondence with the galileon and/or Horndeski Lagrangians of progressively higher order. 
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includes BD. They all correspond to “sub-spaces” of different dimensionality in the space of 
theories, schematically shown in the following figure. 


H45 



Figure 1: Pictorial representation of the space of theories spanned by the EFT coordinates p, P 3 and e 4 . 

The outline of the paper is as follows. We give an overview of the effective held theory 
formulation in Sec. [2j focusing on the quasi-static, sub-Hubble (or Newtonian) limit. We then 
work out the expressions for the gravitational coupling(s) and the gravitational slip parameter. 
We further reduce the degrees of freedom of the formalism by implementing viability conditions 
and by fixing the background expansion history (we choose the ACDM model which best fits 
the Planck’s satellite data !>• In Sec. [3j we introduce the parametrization of the requisite 
time-dependent functions, and systematically scan the parameter space of general dark energy 
theories. We highlight the generic and robust predictions for the relevant large-scale-structure 
observables, such as the linear growth rate or the gravitational slip parameter. We conclude in 
Sec. [4] with a summary of the main results, and qualitative explanations for them. 


2 Effects of dark energy on cosmological perturbations 


There is a range of scales on which extracting perturbation observables from modified gravity 
(MG) theories is relatively straightforward: the window of comoving Fourier modes k s h < 
k < k a \. For momenta less than the non-linear scale, k n \ — (lOMpc) -1 , one can trust linear 
perturbation theory. For momenta well above the sound horizon scale k s h ~ aH/c s ( c s is the 
speed of sound of dark energy fluctuations), one can neglect the time derivatives of the metric 


and scalar fluctuations in the linear equations, the so called quasi-static approximation 56,57 


In the quasi-static regime, it is possible to compute algebraically the effective Newton constant 
G e s(t, k) and the gravitational slip parameter ri(t,k ) of a given MG theory. The entire set of 
perturbation equations then reduces to 


k 2 

- —$ = ATrG eS (t,k)5p ri 
a z 

T 


(2) 

(3) 


where we have adopted the following convention for the perturbed metric in Newtonian gauge, 


ds 2 = — (1 + 2<3 i)dt 2 + a 2 ( 1 — 2 'h)5ijdx l dx^ . (4) 

Eqs. © and © should be supplemented by the equations for the matter fluctuation 5p m . Note 
that we work in the “Jordan frame”—one where all matter fields are minimally coupled to the 
metric (see |27j for a relaxation of this assumption within the effective field theory formalism) 
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and do not couple directly to the dark energy fluctuation. This is the frame of most direct 
physical interpretation 153], where bodies follow geodesic^] 

In the rest of this section, we derive G e g- and r/ in the framework of the effective field theory 
of dark energy (EFT of DE). We also clarify the relation between G e g and the Newton constant 
£rjv measured, for example, in Cavendish experiments. Finally, we summarize the viability 
conditions to be imposed on our parameter space. 


2.1 Gravitational couplings and gravitational slip 


The linear cosmological perturbations of the class of theories considered in this paper are de¬ 
scribed by the unitary gauge action given in App.[Aj By moving to the more practical Newtonian 
gauge, we reintroduce the perturbation it of the scalar field (dark energy fluctuation), along with 
the Newtonian potentials T and T defined in Q. The action quadratic in these quantities is 
very involved, but can be considerably simplified by the following approximations. First, we ap¬ 
ply the quasi-static approximation, and neglect the time derivatives in the gravity-scalar sector. 
Second, in surveys of large scale structure, we generally observe modes well inside the Hubble 
horizon. We thus ignore mass terms in the perturbation quantities it, $ and \k, because they 
are naturally of order Hubble and are small compared to the gradient terms—we retain the 
lowest gradients in the spirit of a gradient expansion)^} In this limit, all the dark energy models 
containing up to one scalar degree of freedom, and whose equations of motion have no more 
than two derivatives, are described by [35] 


5= / aM 2 


(VT) 2 - 2(1 + e 4 )V$Vf - 2(n + e 4 )VTV7r 


+ {p - ^ 3 )V$Vtt - ( C + y - He± + Hh ) (Vvr) 2 


(5) 


a Qdpm, 


where 5p m is the perturbation of the non-relativistic energy density, a dot means derivative 
w.r.t. proper time and n represents the perturbation of the scalar field. The non-minimal 
coupling p and the “bare Planck mass” M are related by 

dlnM 2 (t) 

■ <6) 

and the expression of C is given as 

C = (T) 

In the above, p m oc a^ 3 4 is the background energy density of non-relativistic matter. In minimally 
coupled scalar field models, we can think of C as related to the kinetic energy density of the 
field, C ~ 0 2 /M 2 . To simplify the notation, we have also defined with a circle some “generalized 
time derivatives”: 


p 3 = p . 3 + pp 3 + Hp 3 , 
£4 = e 4 + pc 4 -j- H 64 . 


( 8 ) 

(9) 


Variation of ([5]) with respect to the curvature potential \k produces the algebraic relation 


'P — (1 + £ 4 )^ + (p + ii)TT 


( 10 ) 


3 Geodesic motion can be violated in certain theories that exhibit screening such as in chameleon theories, 
though this assumption remains valid for the kind of screening found in the galileons 


54 


4 Note that we are not exploring chameleon and f(R) models, which phenomenologically require a mass for 


the scalar field much larger than H 36 (see also footnote 1). 
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that, once substituted back into the action, gives 

S = f aM 2 { — (1 + e 4 ) 2 (V$) 2 + \p - p 3 - 2(/i + e 4 )(l + e 4 )] V$Vvr (11) 

+ (/i + e 4 ) 2 + —— We 4 + He^ (V7r) 2 |' — a 3 <$>5p m . 

We can solve the coupled 7r- < I > system by taking the variation with respect to ^ r, giving 

2 C + 2(/i + e 4 )" + p 3 — 2iW 4 + 2Hl^ vr = [/r — p 3 — 2(/r + e 4 )(l + e 4 )] $ . (12) 

We thus conclude that the Newtonian gravitational potential satisfies the Poisson equation 


where 


GW 


87tM 2 (1 + e 4 ) 


k 2 

— ., * 1 1 — ATTG e s(t)p m 6 n 


2 C + p 3 — 2 i/e 4 + 2F/e 4 + 2 (p + e 4 ) 2 

M + h)(p~ Pz) _ (p ~ P3) 2 
1 + e 4 2(1 + e 4 ) 2 


2C + fa — 2 JPe 4 + 2 He^ + 2- 


(13) 


(14) 


As noted, we are neglecting any possible scale-dependence of G e g. The fc-dependent correc¬ 


tions to (14) become important at large distances, how large depending on the size of the mass 
terms that we have neglected in (40). If the scalar degree of freedom, as we are assuming, plays 


a relevant role in the acceleration of the Universe, its mass is expected to be of order Hq or 
lighter. Only Fourier modes approaching the Hubble scale are affected by such mass terms. 

We need now to specify the relation between the effective Newton constant GW and the 
standard Newton constant G]\r —or, equivalently, the Planck mass Mp\. Powerful Solar Sys¬ 
tem |60,61 and astrophysical |62,631 tests impose stringent limits on modified gravity. Re¬ 
alistic models must incorporate “screening” mechanisms that ensure convergence to General 
Relativity on small scales and/or high-density environments. Horndeski theories mainly rely 
on the Vainshtein mechanism, which is now quite well understood also in a cosmological time 
dependent setup [48-51 . In the vicinity of a massive body, the gravitational contribution of the 


scalar field fluctuation it is suppressed because non-linear scalar self-interaction terms locally 
change the normalisation of the 7 f 2 term in the quadratic action 45 . This is equivalent to 


switching off the couplings of the canonically normalized 7r to the other fields T and <f>. Such 
an astrophysical effect cannot be encoded in our quadratic action ([5 ), whose couplings depend 
only on time and not on spac^J However, by inspection of action (Tl]) we conclude that in a 
screened environment, the “bare” gravitational coupling M~ 2 gets simply dressed by a factor 
of (1 + e 4 ) -2 . This is what we obtain in the Poisson equation if we switch off the mixing term 
V4>Vvr. We thus define a screened gravitational coupling 


Gsc(t) = 


8 vrM 2 (l + e 4 ) 2 


valid only in (totally) screened environments. 

Since we live and perform experiments in a screened environment, the value of G c 
today is the Newton constant measured for example by Cavendish experiments, 

1 1 


(15) 


evaluated 


G n = 


87rMpj 


87 rAf 2 (f 0 )[l + e 4 (t 0 )] 2 


(16) 


5 There’s also a “temporal” or “cosmological” screening which screens out modified gravity effects at high 
redshifts 38 : even though our EFT formulation does not provide a microscopic description of how this comes 


about, it does effectively account for such a temporal effect via the time-dependent functions, 71 , P 3 and £ 4 . 
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To complete the discussion, we would like to mention that the coupling of gravitational 
waves to matter contains, instead, only one factor of (1 + € 4 ) at the denominator, and not two 
as G sc 21,35 . In ref. 37], the corresponding mass scale has been associated with the “Planck 
mass”. On the basis of the above arguments, here we find it more natural to define Mp\ as 
in (16) instead. 

In summary, beyond the “bare” mass M multiplying the Einstein Hilbert term in the unitary 
gauge action (40), we can define three gravitational couplings: 


G gw = g 7 r TVi' 2 ( 1 -t-e 4 ) * s cou Pli n g °f gravity waves to matter. 


G e g given in eq. (14) is the gravitational coupling of two objects in the quasi-static 
approximation and in the linear regime, i.e. when screening is not effective. This is 
the quantity which is relevant on large (linear) cosmological scales. 


• G sc = i s the gravitational coupling of two objects in the quasi-static 

approximation when screening is effective. Since the solar system is a screened envi¬ 
ronment, this is also the Newton constant measured by a Cavendish experiment, once 
evaluated at the present day: Gjy = G sc (to). 


Finally we note that even in the absence of anisotropic stress, the scalar and gravitational 
perturbations and T are not anymore of equal amplitude, as predicted by general relativity. 
It is customary to describe deviations from the standard scenario by defining the gravitational 


slip parameter 77 = |), which can be easily expressed within our formalism. By using (10) 


and (13) we obtain 


77 = 1 — 


(/i + £ 4 )(/r + 7/3 + 2 I 4 ) — 64 ( 2 C + 773 — 2He^ + 2Hef) 
2C + fi 3 — 2He± + 2 HI 4 + 2(/r + I 4) 2 


(17) 


As for G e ff, we note that not all deviations of 77 from 1 are screenable in the conventional sense. 


From eq. (10), we see that in the presence of a non-vanishing £ 4 , the Newtonian potentials T 


and <h are detuned from each other, even in environments where the vr fluctuations are heavily 
suppressed, ry sc ~ 1 + £ 4 . 

The above dark energy observables, G e g and 77 , depend on the six functions H(t), M(t), C(t ), 
773 (f) and £4 (t), which are constrained by Eqs. ©> and 0 . Note also that the coupling 
772 (f) does not appear explicitly in the observables, it only plays a role in the stability conditions 
and the speed of sound of dark energy, as we show in the following. Below, we discuss how 
to reduce the dimensionality of this functional space to three constant parameters and three 
functions. 


2.2 Theoretical constraints: viability conditions 

One would expect a theory such as the EFT of DE, which depends on five free functions of time, 
to be virtually unconstrained and therefore un-predictive. What we will show instead is that, 
one can indeed bound the time evolution history of the relevant DE quantities such as G e fr and 
77 defined above. The key is to demand that the DE theories be free of physical pathologies. 

We demand that a healthy DE theory satisfies the following four conditions: it must not 
be affected by ghosts, or by gradient instabilities; the scalar as well as tensor perturbations 
must propagate at luminal or subluminal speeds 47 . In what follows we simply refer to all 
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these criteria as the stability or viability conditions. These conditions must be satisfied at any 
time t in the past, while we do not enforce their future validity. We collect here the algebraic 
relations that enforce these requirements, allowing us to bound the time dependent couplings. 
An in-depth derivation and discussion of these relations can be found in paper I. 


B 

A 


4 = 


where we have defined 


1 + 64 


> 0 

no-ghost condition, 

(18) 

> 0 

gradient-stability condition, 

(19) 

< 1 

scalar subluminarity condition, 

(20) 

< 1 

tensor subluminarity condition, 

(21) 


A — (C + 2^)(1 + £ 4 ) + — P's) 2 , 

B = (C + ^ - hfe 4 + i/e 4 )(l + e 4 ) - 


/i - e 4 


( 22 ) 

(23) 


2.3 Observational constraints: background expansion history 

Recent observations tightly constrain the homogeneous background expansion history of the 
universe - equivalently, its Hubble rate H ( 2 ) as a function of the redshift - to that of a spatially 
flat ACDM model 11,2,41. We thus assume 


(24) 


H 2 (z ) = Hi s 0 (l + zf + (1 - x o) (1 + z) 3(1+w,eff) 


The quantities xo —the present fractional matter density of the background - and w e g —the 
effective equation of state parameter - are free parameters, though observations suggest xo and 
uieff must be close to 0.3 and —1 respectively. Since we are interested in the recent expansion 
history, we have neglected the contribution of radiation. 

The fractional matter density of the background reference model calculated at any epoch, 
x, proves a useful time variable for late-time cosmology, smoothly interpolating between x = 1 , 
deep in the matter dominated era, and its present value xo — 0.3. Its expression as a function 
of the redshift is 

x = - — --tit-• (25) 

x 0 + (l-x 0 )(l + A) 3 ^ff 

Having specified the background geometry H(z), let us now see how this influences the 
functions M 2 (t) and C(t) that are needed to compute G e g (14) and 77 First of all, note 

that M 2 (t) is not an independent free function of the formalism but is related to the non-minimal 
coupling 74 . By inverting eq. ©> we obtain 


M 2 («) = 


Vpi 


(i + «S ) 2 


exp 


dt' p(t'] 


'to 


(26) 


where the initial conditions have been set according to (16). In a similar way, the evolution 


equations for the background (see e.g. [19,24]) result in the expression of C given in 0. There, 
p m represents the physical energy density of non-relativistic matter. By “physical” we mean the 
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quantity appearing in the energy momentum tensor. It scales as a 3 since we are in the Jordan 
frame and p m — 0. With this quantity, we can define the physical fractional energy density 


r>0 _ Pm{t o) 

m - 3M^H% ' 


(27) 


In principle, one could try to measure p m (to) by directly weighing the total amount of baryons 
and dark matter, for instance within a Hubble volume. It is worth emphasizing that, in theories 
of modified gravity, needs not be the sam^] as xq. The latter is a purely geometrical 
quantity, a proxy for the behavior of H(z). It proves useful to define, as in paper l[(J 

k = (28) 

Xq 

In conclusion, the behavior of G e g and p is completely specified once the three parameters 
xo, w e g, and n. and the three functions, p(t), P 3 (t), £ 4 (t), are supplied. Effectively, these 
functions are coordinates in the parameter space of modified gravity theories. 


3 General predictions 

In this section we describe the redshift scaling of interesting cosmological observables such as 
the effective Newton constant G e ff, the linear growth rate of large scale structure / and the 
gravitational slip parameter p. 


3.1 Methodology 

Before presenting our results, let us summarize the framework within which our predictions 
are derived. First of all, we enforce that viable modified gravity models reproduce the ACDM 
background history of the universe. In other words, we choose as background H(z) that of a 
flat ACDM universe with parameters set by Planck (xq = 0.314 and = — 1). 

As for the perturbed sector of the theory, we make the following assumptions 


We adopt a; as a useful time variable interpolating between the matter dominated era 
(x = 1) and today (x = xq ~ 0.3), and expand the coupling functions up to second order 
as follows: 


p(x) = H(l-x) + p^\x — Xq) + pf\x - Xq ) 2 
pz{x) = H(l-x) pf^+p^(x-xo)+pf\x-xo ) 2 
e 4 (x) = (1 - x) pf ] + p^ (x - X 0 ) + p ( 4 2) (x - X 0 ) 2 


(29) 

(30) 

(31) 


On the other hand, we set [12 = 0 since this function’s only effect is that of lowering the 


speed of sound of dark energy (see eq. 20 ) and thus reducing the range of validity of the 


f> One might argue that, since the background evolution is anyway degenerate between the two dark compo¬ 
nents, the distinction between xq and is merely academic. It is true that, in the Friedmann equations, we 
can always relabel some fraction of dark matter as “dark energy”, but at the price of changing the equation of 
state of the latter. By assuming that w e fT be constant—and, in the rest of the paper, u> e ff = — 1 in particular—we 
break the dark-sector degeneracy. Therefore, the distinction between xo and flj), is in principle important. 

7 Note the change of notation with respect to Ref. 24 . defined here is not the present value of the 


function P m defined there. On the other hand, k represents the same quantity. 
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quasi static limit, an approximation used in the present analysis. Secondly, we want to 
recover standard cosmology at very early times. This is why the Taylor expansions are 
multiplied by an overall factor of (1 — x) which serves to switch off the couplings at early 
epochs. This ties the modification of gravity to the (recent) phenomenon of dark energy. 


In the same spirit, we demand that the function C, appearing in G e ff and 77 , grows at early 
epochs less rapidly than H 2 i.e. C/H 2 —> 0 at high redshift. In the limit of an (uncoupled) 
quintessence held (j>, C is nothing else than the kinetic fraction of the energy density of 
the held, C ~ cp 2 /Mpp We are thus requiring that dark energy be subdominant at high 

By inspection of equation 0 ’ 


redshift, which is usually the case in explicit models 38 


we see that we need a cancellation between the last two terms on the RHS. By substituting 
p m = SMp^KxH 2 and using ( [28] ), we see that this condition is equivalent to setting the 
value of M 2 at early times J^] 

M 2 —> . (32) 


By using (26), this imposes a constraint on the Taylor coefficients of p: 


In [rc(l + e°) 2 ] = 


1 - so 

6 w e ?£ L 


2p[ 1} +Pi ] (1 - 3x 0 ) 


hr (x 0 ) 
3w eS 


' ( 0 ) ( 1 ) , 2 ( 2 ) 
p\ - S0P1 + X 0 p\ 


• (33) 


As we point out in the rest of the paper, some of the general features that we find are 
indeed related with this requirement. 


• All of our analysis is carried out by assuming k = 1, i.e., by requiring that the physical 
and the effective reduced matter densities coincide at the present time. This choice is 
somewhat suggested by observations. Direct observations of the mass-to-light ratio on 
large scales lead to a value of f(see e.g. 139]) that is compatible with the “geometrical 
quantity” (here called xq) obtained by fitting cosmological distances. 


We explore the space of non-minimal couplings p, p% and £4 that covers the entire set of 
Hordenski theories by randomly generating points in the 8 -dimensional parameter space 
of coefficients p$ (note that in the p -sector the three are related by eq. 33). We reject 
those points that do not pass the stability conditions (18)-(21) until we have produced 10 4 
viable models. Note that, when dealing with the full 8 -dimensional volume of parameters, 
the chance of hitting a stable theory is much lower than 1%. Because we want the quasi¬ 
static approximation to have a large enough range of applicability, we also impose c s > 0 . 1 , 
which allows us to cover the Fourier volume of the EUCLID mission [57j. We emphasize, 
however, that this is a very weak selection effect on our randomly generated models. 


• There is no “natural unit” in the space of Taylor coefficients. Mostly, we randomly 
generate them in the interval p$ E [—2,2], so that the relevant observables such as 
fa 8 (see Figure 4) span the interval of uncertainty of current data. We have checked the 
robustness of our findings by expanding the ranges for different coefficients, and by raising 
the order of the expansions. A sample of these checks is reproduced in Figure[3} Additional 
checks that we have done include modeling the EFT couplings as step functions of variable 
height and width at characteristic redshifts and using different parameterizations for the 
BD sector. In particular, the one used in paper I is illustrated in App. [B] 

s This is the same relation found in paper I with a slightly different reasoning. 


10 















3.2 Evolution of the effective Newton constant 


As mentioned, the redshift dependence of the effective Newton constant appears to be rather 
constrained by the stability conditions. Indeed, Figure [2] shows that the effective Newton 
constants in our large class of theories have similar evolution pattern over time. The universal 
behavior of G e g(x) is best captured by the following few features. First, all the curves display 
a negative derivative at x = 1, which implies stronger gravitational attraction (G e g > Gn) at 
early epochs (z > 2). This behavior was proved analytically in paper I for a restricted set of 


models in which the Taylor expansion (29)-(31) is retained only at zeroth order. By inspection, 


we see here that the effect is still there at any order in the expansion. The amplitude of such an 
initial bump in G e g(x) varies, with a few models displaying also conspicuous departures from 
the standard model. Also the width of the time interval over which this early stronger gravity 
period extends is quite model dependent. We find that the bump can be somewhat leveled out 


by giving up the requirement that DE be subdominant at high redshift (eqs. 32 and 33). 

Even more interesting is the fact that all the models consistently predict that the amplitude 
of the effective gravitational coupling is suppressed G e g < Gn in the redshift range 0.5 < 
z < 1 before turning stronger, once again, at around the present epoch. The reason for the 
intermediate range suppression is that at those redshifts the dominant contribution to the total 
of G e g is given by the screened gravitational coupling G sc . Stability conditions always make 
G sc lower than Gn during the whole evolution. The characteristic S-shape pattern shown in 
Figure [2] is common to all models and does not depend on the degree of the Taylor expansion 
adopted for the coupling functions (29)-(31). 


We should mention that, within covariant Galileon theories—a subclass of the models consid¬ 
ered here—the same qualitative behavior of G e ff was found (see e.g. Ref. 158j, Fig. 9, Ref. [59] , 
Fig. 3), although the background evolution in that more constrained case is different from 
ACDM. Regarding the weaker gravitational attraction at intermediate redshifts, our results are 
in agreement with those of the recent paper [641. Less clear to us is the role played by the 
gravitational wave speed ct (related to e 4 , in our language) - the author of |64| claims ct < 1 is 
important for obtaining a weaker gravitational attraction. We find that a weaker gravitational 
attraction can be—and, in most cases, is—achieved also in theories ( BD , H3 ) where ct = 1 
(see also j;65j). 

A way to make sense of why the effective gravitational constant is stronger/weaker than 
the corresponding standard model value at characteristic cosmic epochs, is to decompose its 
amplitude into two distinctive parts. We can think of G e g/G]\r as the product of two terms, 


G e ff 

G n 


G e g G sc 
Gsc Gn 1 


which can be expressed as (see eqs. 14 and 22) 


G e ff 

~GZ. 


= 1 + 


1 + e 4 (n- 11 3 


B 


1 + 64 


— (/r + e 4 ) 


Mpi 


G s 


G N M 2 (t)(l + e 4 ) 2 


(34) 


(35) 


respectively. Stability conditions (19) and (21) imply G e g/G sc ^ 1. Physically, this means that 
the scalar field contribution to the gravitational interaction is always attractive, as expected 
from a (healthy) spin-0 field. This circumstance is displayed in the second column of Figure [2j 
The behavior of G sc /Gn , the last column of Fig. [2j also has a physical interpretation related 
to the stability of the models, although somewhat more subtle. First, note that the value of 
such a quantity today is unity by definition—as we have argued in (16), Gn = G sc (to)— while 


at early epochs (x = 1) it is given by condition (32). On the other hand, the overall behavior of 
Gsc/Gtv as a function of time can be understood as a product of the two independent factors, 
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Figure 2: The behavior of G e ff (left panels) and its two separate contributions, from eq. (34) (central 
and right panels), as a function of the reduced matter density x is shown for a sample of 10 a randomly 
generated models of viable BD, H3 and UJ5 theories. The dotted vertical lines identify, from left to 
right, the cosmic epochs 2 = 0.5, z = 1 and z = 2. The thick red line represents ACDM. Note that H3 
and H4-5 do not seem to approach this limit in a continuous way. The point is that BD , H3 and H^5 
correspond to “subspaces” of progressively higher dimensions in the theory space (see figure in Sec. [Tj). 
ACDM represents a well-defined limit for all models, but the volume of stable theories for H3 and HJ^5 
asymptotically reduces to an “hyperplane” of lower dimensions in the vicinity of ACDM. As a result, the 
neighborhood of ACDM is effectively a volume-zero subset for our random generator of theories. The 
Monte Carlo procedure does recover ACDM as a limiting case of BD nonetheless. 


M~ 2 and (l + e 4 )~ 2 . The latter quantity is always lower than unity because tensor perturbations 
are assumed to propagate at subliminal speed (see eq. 21). Also M~ 2 decreases as a function 
of the redshift (i.e. backward in time) at around the present epoch. The physical reason is 
better understood in the “Einstein frame”—the frame in which the metric is decoupled^ from 
the scalar—which for the background evolution simply reads gffi r^j M 2 (t)g^ (see e.g. 1191). A 


growth of M as a function of the redshift means less acceleration in the Einstein-frame, thus 
implying that the observed acceleration in the physical (Jordan) frame is due to a genuine 
modified gravity effect (self-acceleration). Therefore, the third column of Figure [2] provides a 
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Figure 3: Percentage of randomly generated EFT models with growth rate / lower than that predicted by 
ACDM Planck-cosmology as a function of the redshift 2 . We have checked the robustness of our analysis 
by changing the interval of the randomly generated Taylor coefficients (center) and by augmenting the 
order of the Taylor expansion (29)-(31) (right). 


rough estimate of the amount of self-acceleration for the various randomly generated models. 
Curves that deviate the most from ACDM (the red straight line) represent models with strong 
self-acceleration, while the opposite cases represent models in which acceleration is essentially 
due to a negative pressure component in the energy momentum tensor. Our EFT approach 
covers both types of behavior in a continuous way, although it will be interesting to understand 
the specific features of those models that are truly self-accelerating 
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3.3 The growth rate of matter fluctuations 

The universal evolution of the effective Newton constant is expected to result in a characteristic 
growth history for the linear density fluctuations of matter. The prediction for this observable 
is obtained by solving the equation 


3w e ff(l - x)xf'(x) + f(x) 2 + 


2 


3 

2 


(Veff (1 - X) + 1 ) 


/(*) 


_ 3x G e ff 


(36) 


which approximates the true evolution in the Newtonian regime (below the Hubble scale) and 
well after the initial, radiation dominated phase of cosmic expansion. Here, / = din 5 m /din a, 
where 5 m is the matter fluctuation and a is the scale factor. 

Figure [3] shows that the expectation of a pattern of stronger/weaker growth phases with 
respect to the prediction of the standard model of cosmology is confirmed. Understandably, 
since / obeys an evolution equation sourced effectively by xG e g, it responds to periods of 
stronger/weaker gravity with a time-lag. Moreover, as an integrated effect, / has a smaller 
spread compared to G e ff. 

The first thing worth emphasizing is that essentially all modified gravity models with the 
same expansion history of ACDM consistently predict that cosmic structures grow at a stronger 
pace, compared to ACDM, at all redshifts greater than 2 ~ 2. 

A second distinctive feature is that non-standard models of gravity are generally less effective 
in amplifying matter fluctuations during the intermediate epochs in which cosmic acceleration 
is observed, i.e. in the redshift range 0.5 < 2 < 1. Figure [3] shows that 95% of the growth rates 
predicted in the BD , H3 and H^5 classes of theories are weaker than expected in the standard 
ACDM scenario. 
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Figure 4: The redshift evolution of fas expected in 10 3 stable EFT models is shown and compared to 
data. Error bars represent 1 <t s.d.. The evolution predicted by the best Planck fit is shown in red. 


These predictions can be compared with observations. A collection of available measure¬ 
ments of the growth related quantity fag is presented in Figure [4] and compared to the most 
generic predictions of BD, H3 and Hf5 theoriceJ^J The current errorbars are still large, thus one 
should not read too much into these plots. Nonetheless, it is intriguing that the data suggest less 
growth than is predicted by ACDM. If this holds up in future surveys, it would be important to 
check that growth is stronger than ACDM at 2 > 1, as is predicted by the bulk of the models. 


3.4 The gravitational slip parameter 


While the peculiar velocity of galaxies falling into the large scale overdensities of matter con¬ 
strains the possible growth histories of cosmic structures, CMB and weak gravitational lensing 
provide complementary probe of gravity, notably they allow one to test whether the metric 
potentials and T are indeed equal, as predicted by GR in the absence of anisotropic stress, 
or differ as predicted by most non-standard models of gravity. 

Ref. [7] carried out an interesting joint analysis of current data using these three probes, 
obtaining constraints on two characteristic observables that are sensitive to non-standard be¬ 
havior of the metric potentials. These are the gravitational slip parameter rj, already introduced 
in Sec. 2.1 and the light deflection parameter or lensing potential X, defined in Fourier space 
by 

k 2 

-2 ('F + $) = 87rX(t, k)p m Sm , (37) 


which can be expressed as 


X = 


GW 1 + 4 
Gn 2 


(38) 


For BD -like theories it is straightforward to show analytically that the amplitude of the 
curvature potential T is never greater than the Newtonian potential <h, that is, at any epoch, 
rj(t) < 1. For this specific class of theories, the lensing potential reduces to X = G sc /Gpf. From 


9 We rescale the Planck value of the present day rms amplitude of matter fluctuations on a scale 8/i _1 Mpc 
as ag FT (z) = -j^-^xuTTM ag CDM (z = 0) where D+ is the growing mode of linear matter density perturbations 
in various dark energy models, and where the present day value ag CDM ( 0) is set to the Planck value 0.834. 
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Figure 5: The amplitude of the gravitational slip 77 and lensing potential E as a function of the reduced 
matter density x for 10 3 EFT models. The expected values of both these functions in a ACDM cosmology 
are shown in red. The vertical dotted lines correspond from left to right to z = 0.5, z = 1 and z = 2. 


our earlier results (see Figure [2j right panels), we see that this observable cannot be larger than 
unity at any cosmic epoch, and must be equal to 1 (the ACDM value) at the present time. 

When additional degrees of freedom are allowed (like in H3 and/or H/5 models), we still 
find distintive features in the evolution of 77 and E. Indeed, the slip parameter is always smaller 
than unity at any redshift except possibly in the window 0.5 < z < 1, where relevant deviations 
from the ACDM expectations can be observed. Moreover, 77 is never larger than ~ 1.5 at any 
cosmic epoch. Similar to the case of HD-like theories, the lensing potential S is weaker than the 
standard model value at high redshifts (z > 0.5-1), but becomes stronger (greater then unity) 
in recent epochs. Indeed, virtually all H3 and H45 models predict an amplitude of E greater 
than 1 at the present time. 


3.5 Comparison with ad-hoc phenomenological parameterizations 

Given the lack of a compelling model to rival GR, it is common to parameterize deviations 
from GR in a phenomenological manner j7 31,68 . Previous works used simple monotonic time 
evolution for the relevant large-scale structure observable 0(t) (be it G e g/G]y, E or 77 ). For 
instance, the time evolution is often chosen to be proportional to the effective dark energy 
density implied by the background dynamics, 


0(f) = 1 + O 0 (l - x). 


(39) 


so that the modified gravity scenario converges to the standard picture at high redshifts. Note 
that such a simple ansatz corresponds to straight lines in our figures [2] and [5j a time dependence 
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Figure 6 : Left and central panels: Correlations between G e g/G n — 1 and 77 — 1 (top line), Geff/Gjv — 1 
and E — 1 (bottom line) at redshifts 1 and 2 for 10 4 stable Hf5 models. Right panel: The locus of points 
of the stable modified gravity theories encompassed by the EFT formalism is also shown at redshift equal 
to 0. The 2D likelihood constraints obtained by Planck |7 using a phenomenological parameterisation 
for the time evolution of E and 77 similar to that shown in eq. (39) is overplotted. 


that does not seem to correspond to any viable single scalar field model. Instead, it appears 
from our analysis as if there is an additional time scale, roughly located between z = 0.5 and 
z = 1. All relevant quantities undergo a transition at that scale, typically reaching a local 
maximum ( 77 ) or a minimum (G e g/Gi\r and E). This behavior is not captured by the “flat” 
parameterization (39). 

More interestingly, once the EFT formalism is fully and consistently applied, G e g/G^ and 
77 —as well as G e g/G^ and E—appear not to be functionally independent quantities. A strong 
correlation exists among these pairs of observables which, in addition, evolves as a function 
of time. This correlation is displayed in Figure [ 6 ] at three different redshifts for the general 
class of H45 models. This is yet another instructive example of the general property already 
discussed in paper I, i.e. the requirement of physical stability significantly narrows the region of 
the parameter space. Despite the large number of degrees of freedom encompassed by the EFT 
formalism, model predictions are tightly confined within a limited region in the 2 -dimensional 
planes G e g/G]\r — 77 and G e g/GN — E at a given redshift. 

Existing constraints on G e g/G]\r, E and 77 should be interpreted with care: they are sensitive 
to assumptions made concerning their time evolution. This is illustrated in the last panel of 
Figure [bj We superimpose the observational constraints on these variables set by Planck ['Tj— 
obtained by assuming the ansatz (39) for their time dependence—with the final (at z = 0) 
values of the same quantities simulated with 10 4 randomly generated viable theories. Fig. [ 6 j 
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interpreted at face-value, suggests viable EFT models—with their characteristic patterns shown 
in Figures [2] and [5]- span a region of the parameter space that is mostly excluded by data. Such 
an interpretation is unwarranted. Suitable (non-monotonic) time-dependence of G e g/Gi\r, £ and 
rj, in accordance with predictions from viable EFT models, must be included when performing 
the likelihood analysis [701. 


4 Discussion and conclusions 


We explore the phenomenological implications of a large class of dark energy models that result 
from adding a single scalar degree of freedom to the Einstein equations with equations of 
motion at most second order and allowing for non-minimal coupling with the metric. To do so, 
we exploit the effective field theory (EFT) of dark energy. It describes with a unifying language 
different, apparently unrelated, modified gravity theories, and allows us to directly connect these 
theories to observational signatures of departure from general relativity. By enforcing stability 
and (sub-)luminality constraints, we ensure whatever observational signatures that emerge are 
from theories free of physical pathologies. 

In paper I the effective field theory of dark energy was used to predict the linear growth 
index of matter density perturbations. In this paper, we extend the formalism to compute a 
larger set of perturbation observables, including the effective Newton constant, the linear growth 
rate of density fluctuations, the lensing potential and the gravitational slip parameter. They 
are expressed in terms of the structural functions of the EFT formalism. Physical viability 
(i.e. freedom from pathologies) places strong constraints on these structural functions, greatly 
enhancing the predictive power of this approach. 

In particular, the time dependence of G e g(f) (the effective gravitational coupling in the 
Poisson equation 13), /(f) (the growth function; eq. [36]) , and 77 (f) (the slip parameter; eq. © 


for all Horndenski theories can be economically described in terms of three constant parameters 
(xq, w e fi, and k: eqs. 25, 28) which control the evolution of the background metric, as well 


as three functions (yu(i) , /u .3 (t) and 64 (f); eq. [5]) which, being active in the perturbation sector, 
determine how matter fluctuations evolve in time at linear level. In this study we are interested 
in modified gravity models that are kinematically equivalent, i.e. share the same expansion 
history, but differ in their dynamical properties, i.e. predict different growth histories for large- 
scale structures. We thus factor out the background contributions to the amplitude of the 
large-scale structure observables by simply fixing the three constant parameters xo,w e s, and k. 
We do this by requiring that the background expansion rate H(t ) be that of a ACDM model, 
notably the model that best fits Planck data [2|. The space of modified gravity theories is thus 
generated by the three non-minimal coupling functions 77 (f), A* 3 (i) and 64 (f). To proceed further, 
we Taylor expand each EFT function, and scan the theory space by generating the expansion 
coefficients in a Monte Carlo fashion. An important point is that not all the models produced 
in this way are viable. There is a set of stability conditions that must be satisfied in order 
for a modified gravity theory to be healthy. These act as a severe selection on the parameter 
space and allow us to identify distinctive qualitative features of single scalar field dark energy 
models. We consider 10 4 models surviving the stability criteria while matching the same cosmic 
expansion history. 

This paper complements paper I in providing a maximal coverage of the viable theory space, 
and strengthening the robustness and generality of the predictions. We confirm two central 
findings of paper I, i.e. a) the space of theories naively allowed by cosmological data is much 
reduced once the viability criteria are applied; b) the vast majority of stable Horndeski theories 
produce an overall growth of structure at low redshifts that is weaker than that expected in a 
ACDM model. Here is a summary of our most important findings. 
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• Modified gravity does not necessarily imply stronger gravitational attraction. There are 
two effects that work in opposite directions. One effect is best seen in the Einstein frame— 
the frame where the metric is de-mixed from the scalar—the addition of an attractive 
scalar force indeed makes gravity stronger. This is equivalent to the statement that 
G e s(t) > G sc (t ) (see middle panels of Fig. [ 2 ]). Note that this involves a ratio of the 
effective G (that includes the effect of the scalar) and the “Einstein” G (i.e. G sc ) at the 
same time. A separate effect is encoded in the time-evolution of the bare Planck mass 
M(t ) (or of Gsc). Models that achieve cosmic acceleration by virtue of a genuine modified 
gravity effect do so by giving G sc a smaller value at intermediate redshifts (z ~ 0.7) 
than the current one, i.e. G sc (t ) < G sc (£o)- When comparing structure growth in such 
models against that in general relativity, we are really asking the question: what is the 
ratio G e ff(i)/G sc (io), since G sc (£o) is the value of the gravitational coupling observed 
in the solar system today (we call this Gjy). The ratio G e ff(t)/G sc (to) can be written 
as [G e fr(£)/G sc (£)] x [G sc (t)/G sc (£o)]> hence it is a competition between two factors one 
larger, and one smaller, than unity. Indeed, we find that in a narrow redshift range around 
z — 0.7 virtually all MG models display G e fr < Gn = G sc (to)- This feature is present also 
in the growth function /, though to a lesser degree because / is an integrated quantity, 
dependent also on the behavior of G e s at earlier cosmic times. In the redshift range 
0-6 < z < 0.8, only l%(/5%) of modified gravity models predict a value of /(/ fas) that is 
larger than the ACDM value. Recent observations suggest that the growth of structure in 
the low redshift universe is somewhat weaker than is expected in the best-fit ACDM model. 
It remains to be seen if this is confirmed by higher precision measurements in the future. 
It is nonetheless interesting that a weaker growth is predicted, at intermediate redshifts 
around 0.7 or so, by the EFT of dark energy without fine tuning of the parameters. 

• We identify an epoch of super-growth (i.e. faster growth than in ACDM) for cosmic struc¬ 
tures. At sufficiently high redshifts, it is the first factor in [G e ff (£)/G sc (£)] x [G sc (t)/G sc (to)] 
that wins over the second factor. At high redshifts, there is no need for cosmic acceleration, 
thus G sc (t) /G sc (fo) returns to unity, while G e s(t)/G sc (t) remains larger than unity due to 
the presence of the scalar force. Indeed, as already demonstrated analytically in [24| for a 
more restricted number of models than is considered here, dG e ff(z)/dz is negative at the 
big-bang, implying a stronger gravitational attraction at high redshifts. More concretely, 
we find that for redshifts greater than z = 2, both / and fas are always larger than their 
ACDM values. This prediction could be used to potentially rule out the whole class of 
EFT models investigated in this paper. Structure growth at such high redshifts is poorly 
constrained by current data, but the situation will improve with future surveys such as 
Euclid, DESI or SKA. 

• For Brans-Dicke theories, the gravitational slip parameter is at most one: 77 ^ 1. This 

can be shown analytically from Eq. or qualitatively understood as follows. Absent 

anisotropic stress, = T in Einstein frame. The Jordan frame $ and T are obtained from 
their Einstein frame counterparts by a conformal transformation involving the scalar 7 r. 
Because <h and T comes with opposite signs in the metric, the conformal transformation 
affects them differently, leading to r) — 1 = (\k — < f > )/4 ) oc tt/Q, with a proportionality 
constant that makes — 1 ^ 0 (recalling that 7r oc according to the equation of motion; 
see eq. [I2| ). We find that this feature is substantially inherited by all theories, except for 
a narrow redshift range around 0.5 < z < 1 when r/ exceeds one. For all theories we find 
g < 2 at any redshift. 

• The requirements of stability and (sub-)luminality help greatly in narrowing the viable 
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parameter space of dark energy theories. Thus, despite the presence of several free func¬ 
tions in the EFT formalism, the class of theories we study is highly predictive. Moreover, 
the background expansion history is already well-constrained by current data, further 
limiting the form of the free functions. Applying the EFT formalism, we find interesting 
correlations between different large scale structure observables. In particular, at a given 
redshift, there exists only a narrow region in the 2-dimensional planes G e g/G]\r — V and 
GW/G'tv — S where data can be meaningfully interpreted in terms of viable theories. 


The systematic investigation by means of the EFT formalism of what lies beyond the stan¬ 
dard gravity landscape is still in its infancy, and a number of improvements would be desirable. 
For example, it would be interesting to work out the consequences of relaxing some of the 
constraints imposed in our current analysis. How might our predictions be modified if a dif¬ 
ferent initial condition were chosen (k/ 1 ), or if the background expansion history is altered 
(iCeff ~f~ —1)? It would also be interesting to consider early dark energy models, in which a rel¬ 
evant fraction of dark energy density persists at early times (he. not imposing eqs. 32 and 33). 
We have focused on scales much smaller than the Hubble radius in this paper. As data improve 
on ever larger scales, our analysis should be extended to include possible scale dependent effects 
that come from mass terms for n that are of the order of Hubble. Lastly, it would be useful to 
isolate models that have genuine self-acceleration, and compare their predictions with the ones 
studied here (which include both models that self-accelerate and models that do not). 
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EFT action and relation with other equivalent parameteri- 
zations 


While referring to the existing literature 19-23 for more details, here we just display the EFT 

wl 

S = S m [g^ u , ipi] + ld 4 x^g - 


action written in unitary gauge, which implicitly defines the couplings ([!]), 

M 2 {t) 


[R - 2A (t) - 2 C{t)g 


.00 


-l4(t)(Sg 00 ) 2 ~ hs (t)5K5g 00 +efit) 6 K" - 6 K 2 + 


oo 


(40) 


Here the metric is the Jordan one, to which matter fields 'fii minimally couple (we are assuming 
the validity of the weak equivalence principle here, see 27 for a relaxation of this hypothesis). 


The functions A and C are constrained by the Friedman equations (see Sec. 2.3 or Ref. 24). 

In order to arrive at the action ([5j from (40) we need to change gauge. The Stueckelberg 
procedure (see e.g. 123 for details) allows one to exit the unitary gauge by performing a time 
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diffeomorphism t —> t + ir, where n represents the fluctuations of the scalar field, that now can 
appear explicitly in the action and in the equations. The Newtonian gauge then follows from 
further requiring a metric of the form Q. Note that the Stueckelberg procedure does not create 
direct coupling between n and the matter fields, consistent with the fact that we are working 
in the Jordan frame, and the matter action is diff-invariant. 

Finally, it is useful to write a dictionary relating our couplings to the equivalent “a couplings” 
defined in Ref. [37] (see also [26]): 

= M 2 (1 + 64), 

H{ l + e 4 ) H ’ 

2C + 4p 2 

= H 2 { 1 + e 4 ) ’ ( 41 ) 

_ P3 ~ P 
H( 1 + e 4 ) 

_ £4 

1 + £4 


M 2 

OLM 

OLK 

as 

CtT 


B Exploring EFT models by means of a different parameteri¬ 
zation scheme 


In paper I a specific parametrization of the Brans-Dicke non-minimal coupling p(t) was sug¬ 
gested, based on the “physical” equation of state parameter for dark energy w = pd/pd ■ This 
approach provides a complementary way to solve the background/perturbation sector degener¬ 
acy in the action ([ 5 ]) so that the function C is no longer free but is directly linked to p. The 
details about this parameterization scheme (the so called “re-parameterization” in opposition 
to the “/^-parameterization” adopted in the current paper) can be found in Paper I. We briefly 
summarize here only the intermediate steps that are needed to relate these two parameteriza- 
tions, and are thus important for expressing observable such as G e ff and fas in the language of 
paper I. 

A little bit of algebra is enough to show that p and the physical equation of state parameter 
re are related as follows 


p(x) 


■iWeffHjl - X) 
W - rC e ff (1 - x) 


W — W e ff + x{l 


w e g dw 
w dx 


(42) 


One can thus proceed by expanding in Taylor series the function w (instead of directly p, 
as we have done in this paper). 


w(x) = W e ff 


1 - x 0 


1 - kx 0 (1 + el) 2 


+ p^ (x - X 0 ) + p[ 2) (x ~ X 0 ) 2 , 


(43) 


We should mention that the zero-th order term is fixed by requiring that the expansion rate of 
the various EFT models be identical to that of the reference flat ACDM scenario ( u> e fr = — 1 
and xq = 0.314). The other couplings functions are expanded as discussed in section 3.1 of the 
present paper. 

Figure [7] displays some of the results we obtain by implementing the re-parameterization 
scheme. Since there is no natural unit for the coupling functions, we chose to randomly pick 
Pi € [— 0 . 8 , 0 . 8 ] (for w) and p^l £ [— 2 , 2 ] (for p 3 and e 4 ) in order for /<r 8 curves to cover current 
observational data as much as in the ^-parameterization. The resemblance of G e s/G n to the 
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one presented in Figure [2] confirms the robustness of our general findings; the universal behavior 
of G@ff (S-shape) is recovered no matter what parameterization is used for the coupling functions. 
In addition, the fraction of models with lower/higher growth with respect to ACDM displays 
the same general scaling as a function of z. We find that ~ 90% of H45 theories predict 
growth rates below the ACDM prediction in the redshift window 0.5 < z < 1. The change 
of parameterization implies a slight modification of the asymptotic behavior of the coupling 
functions /i and C/H 2 at x = 1, more precisely a modification of the speed at which these two 
functions go to zero during mater domination. This is well illustrated by the tail of the curves, 
the transition to the super-growth epoch is now shifted to a somewhat higher redshift (z ~ 2.5). 
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Figure 7: The plots displayed are computed in the ^-parameterization scheme. The evolution of G e ff(a;) 
(left), fcr s (z) (middle) for a sample of 10 3 randomly generated models of viable H45 theories is shown. 
The percentage of 10 4 randomly generated EFT models with growth rate / lower than that predicted 
by ACDM Planck-cosmology as a function of the redshift z (right) is also shown. The thick red line 
represents the prediction of ACDM. 
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